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^ ! We discuss what the presently collected data tell us about the mechanism of multifrag- 
^ I mentation by comparing the results of two different models, which assume or show an 
^ ' opposite reaction scenario, with the recent high statistics 47r experiments performed by 
QQ ■ the INDRA collaboration. We find that the statistical multifragmentation model and the 
O . dynamical Quantum Molecular Dynamics approach produce almost the same results and 
^ ! agree both quite well with experiment. We discuss which observables may serve to over- 
come this deadlock on the quest for the reaction mechanism. Finally we proof that even 
if the system is in equilibrium, the fluctuation of the temperature due to the smallness of 
the system renders the caloric curve useless for the proof of a first order phase transition. 



1. Models for describing multifragmentation 



^ ' About 15 years ago it has been found that in nucleus-nucleus collisions at intermediate 
energies up to 15 intermediate mass fragments (IMF's) with Z > 3 are created. Since then 
many dedicated experiments have been performed to "nail down" the reaction mechanism 
but despite of all experimental and theoretical efforts we have not found an answer yet. 
The different contributions in these proceedings present a clear evidence for this fact. 
Over the years two different conjectures have been launched and developed, improved 
and further improved without giving a conclusive answer. The reason is, as we will 
see, that the two most advanced but also almost opposite approaches predict the same 
behavior for several key observables. Before we discuss this in detail we will introduce 
shortly the two principal approaches. 

1.1. Multifragmentation as a statistical process 

One may assume that multifragmentation is a statistical process. This means that 
during the reaction at least a subsystem comes to statistical equilibrium and maintains 
this equilibrium during the expansion until it has reached the freeze out volume. Then 
it disintegrates INSTANTANEOUSLY into neutrons, light charged particles (LCP's) and 
IMF's. Because it is assumed that at freeze out the system is in statistical equilibrium 
the disintegration pattern is determined by phase space. This means that one searches for 
all microstates (consisting of fragments and nucleons with a given kinetic and excitation 
energy) which are compatible with a total energy, a total proton and neutron number 
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and a given volume. Each microstate has the same weight and hence one samples all 
microstates in order to obtain jtp^A ■ There are three different models - 0) which 
perform this task. Sequential decay (i.e. the formation of a compound nucleus which 
decays by subsequent fragment emission) has been ruled out because it gives too few 
fragments. 

Unfortunately this task is not as easy as it seems to be, both conceptually and numer- 
ically. Therefore the results of these models differ considerable despite of the fact that 
one is tempted to believe that counting the number of microstates is a well defined task. 

Conceptually there are two problems: 

a) there is no theory which tells us how to treat the unstable states. Usually the number 
of excited states of fragments is calculated via a level density parameter but some of these 
excited states are unstable against particle (usually neutron) emission. Therefore it is 
not evident, whether they should be treated as an excited fragment A or as a fragment 
A-1 plus one neutron. Both descriptions give rise to different microstates. The different 
answers to this problem yield a quite different neutron yield as has been recently found 
out by Toeke and Schroeder 0]. 

b) the freeze out volume can be defined quite differently. Whereas in ref [0 a constant 
freeze out volume is assumed in ref. HI] the freeze out volume depends on the microstate. 
There is no convincing argument why the one assumption should be better than the other. 

Also numerically (especially concerning the Coulomb energy) the approaches differ con- 
siderably. It is the merit of Gross and Sneppen to have discussed all the differences in 
detail and to have demonstrated how it comes that two seemingly identical approaches 
produce quite different results. 

1.2. Multifragmentation as a dynamical process 

On the other hand there are dynamical models. Originally developed to describe proton 
spectra, particle production and the influence of the nuclear equation of state on observ- 
ables like the flow, they have later been further developed to describe multifragmentation 
as well. In 1985 the Quantum Molecular Dynamics model has been introduced which 
describes the time evolution of the n-body Wignerdensity and hence allows for the inves- 
tigation of fragments which are n-body correlations in this context. In this approach the 
nucleons are presented as Gaussian wave functions 



where a characterizes the nucleon. The time evolution of the individual nucleons is 
determined by a variational principle. At the very end a minimum spanning tree is applied 
to identify the clusters. In these models the whole reaction, starting from the initially 
separated projectile and target nuclei until the finally observable fragments and single 
particles is simulated. Besides the Hamiltonian which is, however, well determined by the 
requirement to have stable nuclei with the right binding energy, the only parameters to 
change are the impact parameter and the beam energy. The different final states, created 
in each individual simulation, are caused by the different possibilities to realize the same 
single particle distribution and the same binding energy by choosing different positions 
Xa and momenta Pa for the individual nucleons. For details of this model we refer to the 
references 
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There are two conceptual criticisms of this model. First, that the wave function of the 
nucleons is not antisymmetrized. Therefore one argued, making reference to a Fermi gas, 
that the specific heat will be different as compared to an antisymmetrized system. Hence, 
if the system comes to equilibrium one expected a quite different fragmentation pattern. 
Although this is true of course for a Fermi gas, the Aladin collaboration has found p| 
that this is not true for an interacting system containing nucleons and fragments. The 
above mentioned incertainties of the statistical models have a larger influence on the 
fragmentation pattern than the different specific heats. Second, being a semiclassical 
approach, the fragments have the properties of the Weizsacker mass formula but are not 
real quantum states. In addition effective charges are used. Hence isotopic distributions 
cannot be calculated. 

Nevertheless with two parameters which have to be adjusted, the width of the wave 
function and the fragmentation radius of the minimum spanning tree, the QMD model has 
described many low energy fragmentation data [0,0. At these energies the whole nucleus 
disintegrates into fragments. At higher beam energies, where only the weakly excited 



spectator matter disintegrates into fragments the model in its original version failed 10 



Only recently the reason for this failure has been identified and the model describes now 
spectator fragmentation as well ITU] . 

Careful investigations have shown that in these simulation programs the system passes 



never a state of equilibrium [12|. In almost all cases where multifragmentation has been 
observed the systems is, on the contrary, far away from it. This can easily be seen from 
the videos which have been produced with help if this program ||13||. 



2. Multifragmentation - a theory invariant phenomenon? 

In view of the quite different dynamics assumed or obtained in these both approaches to 
multifragmentation it seems to be astonishing that experiments have not yet decided which 
of the both conjectures is correct. This fact becomes more understandable if one considers 
in similar but somewhat simpler models, what one expects to see. For this purpose we 
assume that the system can be described by canonical thermodynamics, although this is 
not the case, as we will see later. For the question discussed here it serves the purpose. 

First we will suppose a fast fragmentation process as it is the case for a dynamical 
process. If the disintegration of the nucleus is instantaneous each nucleon keeps its mo- 
mentum due to the Fermi motion and one expects an average fragment kinetic energy of 



3/5Ef independent of the fragment size |T^, where Ep is the Fermi energy. The same 



independence one expects if the fragments are formed very late, after the system has been 
expanded while maintaining thermal equilibrium as assumed in the statistical model. This 
requires the opposite, that the disintegration is very slow. Here the average kinetic energy 
is 3/2 T, where T is the temperature at freeze out. Hence it is only the absolute value 
but not the mass dependence of the mean kinetic energy which is different. 

The same ambiguity one finds for the mass yield. Fischer W^ has pointed out that 



finite systems close to a phase transition show a fragment distribution (t{A) oc . But 
this power law dependence appears in many physical processes which are not thermal at 
all |1T6|-[18| . Figure |I| shows schematically this behavior. 



This ambiguity remains valid if one employs the more sophisticated models discussed 
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Figure 1. Similarity of the results one expects to see independent of whether multifrag- 
mentation is a statistical or a dynamical process. 



above. For our studies we use data and simulations for central collisions of the reaction 
50 AMeV Xe + Sn, measured by the INDRA collaboration ||19[. This is the experiment 
where the most precise data are available at the moment. The SMM calculation we are 
discussion about are adjusted to reproduce the subset of events fulfilling the condition of 
completeness (more the 80the total charge and the total momentum are detected) and 
they show a flow angle 9 flow > 60° When comparing with a dynamical calculations, 
we will use the subset of events where the total transverse energy of light charged particles 
Etrans > 4:50MeV. Both models, SMM and QMD, have been extensively used to interpret 
this reaction |^-[^. We would like to mention that an analysis of the QMD calculations 
shows that in this approach the system never passes through a state of thermal equilibrium 
[p^ whereas the application of SMM is only justified if the fragments are formed from an 
equilibrized subsystem. Thus the reaction scenarios in these models are orthogonal. For 
the comparison with the data both calculations have been filtered with the experimental 
acceptance and a centrality cut is employed The increase of the average transverse 



energy of the fragments as a function of their mass, observed in the experiment, is larger 
than predicted in SMM [IS] calculations. Therefore the statistical model calculations 
have been modified by adding a fourth system parameter, a collective energy, which is 
parameterized as Ecou = c * A, where c is a parameter which remains to be determined 
and A is the fragment mass. The best agreement the subset of experimental data fulfilling 
the centrality condition 6 flow > 60° and SMM calculations is obtained with the following 
set of input parameter: 



freeze out density: l/3po 
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source size: 25=78 As = 186 

excitation energy: Ethermai=7-A MeV Ecoh=2.2A MeV 

Even central collisions at intermediate energies have a less or more binary character and 
consequently emission close to beam or target velocity spoils the spectra of particle emis- 
sion from a possible thermal source at rest in the center of mass. Therefore a meaningful 
comparison between statistical model calculations and experimental data is only possible 
around 9cm = 90°. We subdivide therefore the experimental data and the simulations 
into two equal size 27r intervals: 

Bobs : 60° < OcM < 120° 
Bunobs '■ OcM < 60°, 9cM > 120°. 

In Bobs we observe a flat angular distribution and a constant average energy of IMF's and 
LCP's as a function of the emission angle, both being prerequisites for a statistical 



equilibrium. Applying the statistical models, one assumes that in Bunobs a preequilibrium 
component is superimposed to the thermal component. Therefore the statistical model 
cannot be compared to data there. In the dynamical model there is of course nothing like 
a distinction between an equilibrium and a preequilibrium part therefore we can compare 
the results with the data. As fig. ^ shows, QMD describes there the fragment spectra quite 
well |^1[|. Now we come to the comparison in Bobs- In fig- we display on top the average 



kinetic energy of the fragments. On the bottom we show the charge yield distribution. 
We display the results for QMD and SMM calculations in comparison with the INDRA 
data. As one can see, these are well reproduced in both theories, underlining the above 
mentioned observation that these observables are not sensitive to the reaction mechanism. 
Hence we are confronted with the fact that completely different underlying reaction dy- 
namics produces the same key observables for the fragments seen in 60° < 6cm ^ 120°. 
Therefore more complicated observables have to be employed to distinguish between the 
different possible reaction mechanisms. The next more complicated is 



dEdA- 



2.1. Fragment energy spectra 

In fig.^, left hand side, we compare the experimental fragment kinetic energy spectra 
for emission at 60° < 6cm < 120° with the SMM calculation for three different fragments. 
On the right hand side a surprisal analysis is presented to show the structure of the 
deviation between the experimental and theoretical energy distribution. We see that 
the increase of the average fragment kinetic energy as a function of the fragment mass 
- seen is fig. ^ - has a quite different reason in the SMM calculations as compared to 
the experiment. In SMM the maximum of the energy distribution increases due to the 
flow but the high energy slope remains almost constant. This is expected because the 
disintegration due to phase space produces the same slope for all particles. In the INDRA 
data the increase of the average kinetic energy is due to an increase of the high energy 
slope of the spectra. Because the emission of fragments at midrapidity is a rare process 
which is logarithmically suppressed for large charges we do not have the statistics in QMD 
to compare in detail the kinetic energy spectra of fragments with a charge larger than 12. 
For higher charges the fluctuations render the analysis meaningless, for the slopes as well 
as for the spectra. The spectra for selected charges are presented in fig. ^. We see a rough 
agreement of the high energy slope but there are too many fragments at low energy. This 
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Figure 2. Central Collisions (Etrans > 450 MeV). Kinetic energy spectra of fragments 
emitted in forward/backward direction. The energy where fragments have Ebeam/N is 
marked by a line. 



effect is understood but there is presently no remedy to it. It is caused by the nuclear 
interaction range. It is too large as compared to the standard values and screens therefore 
a part of the repulsive Coulomb interaction. Hence one has to conclude that even the 
most complete experiment, available up to now, is not sufficiently precise to distinguish 
in a convincing way between the different theories. The more detailed the comparison is 
made the more the shortcomings of the approaches become apparent. As seen, precise 
high energy fragment data may discard the statistical models because the high energy 
part of the spectra has to have there a mass independent slope. To measured them is, 
however, a very cumbersome task. Therefore it is better to look for a more intelligent 
solution. A possible way to overcome this deadlock we will discuss in the next section. 

In general the statistical model calculations provide excellent agreement with data as far 
as multiplicity distributions, correlations between the fragment masses and other quanti- 
ties, which have nothing to do with the kinetic energy of the fragments, are concerned. 
They fail usually if they are compared with kinetic energy distributions, which show an 
apparent slope temperature of about 15 MeV. Such an high apparent temperature can 
only be created if the excitation energy/N is higher than the binding energy/N and hence 
if the fragments are not stable anymore. Hence at these excitation energies the statistical 
codes produce configurations consisting LCP's only. 
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Figure 3. Mean energies and charge distributions for QMD, INDRA and SMM data in 
60° < OcM < 120° for 50A MeV Xe + Sn, central collisions. 



3. Isospin tracing- a possible method to overcome the present deadlock 

The above discussion has shown that it is all but easy to use the fragment distribu- 
tions, obtained in a single reaction, to determine the reaction mechanism. Therefore it is 
certainly worthwhile to think about a different approach. The idea to use a combination 
of different projectile/target combinations to address this questions, has been advanced 
quite a while ago by the Texas group [^. Changing the isospin of projectile and/or 
target, one may investigate the multiplicity of the different fragments and can compare 
the results with the prediction from a statistical model calculation which assumes that 
the system comes to equilibrium. The problem is, however, that the variation of the 
isospin is rather limited and therefore the effects are tiny. Therefore it is necessary to 
find a convincing analysis of the results. Recently the FOPl collaboration came up with 
an analysis which allows for a direct answer to the question of equilibration Fol- 
lowing earlier work on smaller systems of the Texas group [^] they used an equal mass 



projectile/target combination (Ru and Zr) and performed all 4 possible reactions (Zr-|-Zr, 
Zr+Ru, Ru-|-Zr, Ru-|-Ru). Because the center of mass is the same many detector accep- 
tance problems disappear, and in addition a comparison of Zr+Ru and Ru+Zr allows to 
address directly the question of equilibration as a function of the rapidity. If the system 
comes to equilibrium, in both combinations the same rapidity distribution is expected. 
If, on the contrary, there is a difference in the rapidity distribution, the assumption that 
the particles are emitted from a statistically equilibrated source is excluded. For the first 
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Figure 4. Emission at 60" < 9cm < 120° Energy spectra: we compare the SMM and 
INDRA data. On the right hand side the spectra are displayed, on the left hand side the 
surprisal analysis. 



results of this experiment performed at 400 AMeV we refer to the contribution of Rami. 
Here we only repeat the main result and discuss the interpretation in the framework of 
the models discussed above. In fig.^ we see Rz = "^^z^r^z^n " as a function of the rapidity 
for the 4 different projectile/ tar get combinations. Z is the number of protons observed at 
the given rapidity. Rz gives 1 for Zr+Zr, -1 for Ru+Ru and if equilibrium is reached 
in the mixed combinations. We see clearly that for y^°^ = — ^ — > .25 the protons do not 

ybeam 

show a sign of equilibrium. The same has been observed for the ratio of He^ and triton 
[ P^ . Fig. fright, displays the result of the QMD calculations which have been performed 
by the FOPI-coUaboration. For testing the sensitivity of the results they have employed 
different nucleon nucleons cross sections by rescaling the fitted experimental values by a 
given factor. We see that the proton rapidity distribution depends quite sensitively on 
the cross section. Already an increase by a factor of 2 brings the system close to equi- 
librium. The experimental data are best described by assuming that the cross section in 
the nuclear medium is close to the free one. 

Although this beautiful experiment has shown that at a beam energy of 400 AMeV 
the protons and the light fragments are not emitted by an equilibrated source we have 
no definite answer about the mechanism of multifragmentation yet. Therefore we have to 
continue along this proposition. We have to study the rapidity distribution of fragments 
not only at 400 MeV but at lower energies as well. But this idea to use the isospin degree 
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Figure 5. Emission at 60° < 6cm < 120° Energy spectra: we compare the QMD and 
INDRA data. On the right hand side the spectra are displayed, on the left hand side the 
surprisal analysis. 



of freedom as a tag for equilibration, may free us from the present deadlock and may give 
a final answer concerning the mechanism of multifragmentation. 

4. Canonical models-applicable to nuclear physics? 



In recent times ^^it has been conjectured that the system formed in heavy ion reactions 
can be described by canonical thermodynamics as well. This conjecture is the basis for the 
caloric curve, which displays the temperature as a function of the excitation energy of the 
system. In case that the system suffers a first order phase transition, this function becomes 
flat, i.e. an increase in the excitation energy does not corresponds anymore to an increase 
in temperature because the energy is used as latent heat. The Aladin collaboration has 



claimed to have seen such a sign for a first order phase transition |£3[. The prerequisite 
for the meaningfulness of such a plot is the possibility to measure the temperature and 
the excitation energy at the same time. Since the temperature is extracted from the 
fragment yield it is the temperature of the IMF's which counts. Here we do not discuss 
the question how experimentally these quantities can be measured. We are concerned 
with the question whether for a system as small as that created in nuclear reactions the 
temperature fluctuations for a given excitation energy are sufficiently small to make a 
caloric curve possible. In a large system, the fluctuation of the excitation energy for a 
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Figure 6. left:Rz (see text) as a function of the scaled cm rapidity for the 4 equal mass 
target projectile combinations. Right:QMD predictions for Rz (see text) for different in 
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given temperature is oc where N is the number of particles in the system. In a large 
system, with large N these fluctuations are therefore negligible and for a given excitation 
energy we find a precise value for the temperature of the system. 

In order to discuss this question we would like to start with an explanation of the 
differences between the microcanonical description of a system, as done in the statistical 
multifragmentation models, and canonical thermodynamics. The basis of all statistical 
mechanics is the counting of microstates. To make life easier we take a concrete example. 
Consider an harmonic oscillator whose energy levels i (i > 0) can be occupied by particles 
which have energies = huj ■ i We have suppressed the zero point motion here. The total 
energy of the particles is given hj E = J2 ' ''^i where Ui is the number of particles which 
are on the level i. Now we can calculate the number of different repartition which exists for 
a given number of particles N = J^rii and a given energy E = J2 ^i^i- For example, if we 
have N=2 and E = 2huj we can either have both particles on the level i=l or one particle 
on the level i=0 and the other on i=2. Each of these repartitions is called a microstate. 
g{N, E) is the number of microstates for a given E,N and Q{N, E) = k In g{N, E) is the 
entropy of the system (k is the Boltzmann constant). 

Now, assume that we have two oscillators which can transfer energy and particles and 
have the total energy E = Ei + E2 and the total particle number = A'^i + A'"2. The 
number of microstates of the combined system is given by 

g{N, E)= g2{N,, E,)g,{N -N,,E-E,) (2) 

Ei,Ni 

where gi{Ei, Ni) and (72(^2, E2) are the number of microstates of the oscillators 1 and 2 
for given energies and particle numbers Ei, Ni, E2, N2, respectively. 
The temperature is defined as For our combined oscillators 



g{N, E, E,) = ^ g2{N,, E,)g,{N -N,,E- E,) 



(3) 
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one finds 



kdlng{N,E,Ei) 
dE^ 




(4) 



(5) 
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What can one conclude? Only if g{N, E, Ei) is maximal (at Ei = E^""^) the temperatures 
in the two subsystems are identical. Otherwise the temperatures are different. Hence for 
small systems it is the generic case that both subsystems have different temperatures. 
This is completely counterintuitive because our experience is based on the large systems 
we encounter in real life where - in equilibrium - the two subsystems have the same 
temperature. There the sum in eq. 2 can be very well approximated by the largest term 



Only if this approximation is justified it makes sense to describe the system as a canonical 
system, i.e. in the variables T and /z. Then the system obeys all the laws of equilibrium 
thermodynamics, that the entropy becomes additive and that subsystems in equilibrium 
have the same temperature and the same chemical potential. 

Is this approximation justified in heavy ion collisions where the systems are still rather 
small containing at most around 400 particles? We use the above mentioned micro- 
canonical SMM calculations to find an answer. We identify the two subsystems with 
the fragments with Z > 3 and the LCP's, respectively, and calculate the temperature 
difference between the subsystems. 

In fig. ^ top, we display the number of events (left) and the average number of IMF's 
(right) as a function of the total number of nucleons entrained in the fragments and the 
total fragment energy. We see a rather broad distribution. In the middle we present the 
difference of the temperature and the chemical potential in the two subsystems (LPC's 
and IMF's)as a function of the total mass and the total energy of all IMF's. For calculating 
the temperature difference we have assumed that ^ — ^ t^+'^t^/a "^^ere T is assumed 
to be 7 MeV. We see a rather broad distribution of the chemical potentials and of the 
temperatures. The probability distribution of fragments to be emitted in a microstate 
which shows a | AT | and a chemical potential difference of | A/x |, respectively, between 
the two subsystems is plotted in the bottom row. The chemical potential | AyU | in 
units of the temperature fluctuates by about to 40% around the mean value whereas 
the variation of the distribution of the temperature difference between the subsystems is 
about 1.7 MeV. More precisely the temperature in the subsystem of the fragments cannot 
be determined. Consequently, the temperature values in the caloric curve have to get 
an error bar of this size which renders any conclusion about a phase transition useless 
because one cannot distinguish whether T increases slowly as a function of the excitation 
energy (cross over) or whether it stays constant (first order phase transition). 

This variance of the temperature in the subsystems renders as well the application of 
canonical thermodynamics to heavy ion reaction useless. There is nothing like a well 



(?2(iVi, Ei)gi{N -Ni,E- E^) ^ g^iN^, ErigiiN -Ni,E- Eri 
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Figure 7. Top: Number of events and number of fragments as a function of the total mass 
of all fragments and of their total energy. Middle: Distribution of the temperature dif- 
ference and the difference of the chemical potential. Bottom: Distribution of the absolute 
value of the temperature difference distribution and that of the distribution of the differ- 
ence of the chemical potentials between the subsystems formed by fragments and LPC's, 
respectively, for a SMM calculation adjusted to describe the INDRA results 
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